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This  Semi-Annual  Technical  Report  provides  an  update  on  progress  in  the  analysis 
of  the  seismic  data  collected  during  the  "NGENDEI"  Experiment  sponsored  by  DARPA  in 
early  1983  in  the  southwest  Pacific.  The  experiment,  which  took  place  on  the  Deep 
Sea  Drilling  Project  Leg  91,  was  designed  to  test  the, Marine  Seismic  System  (MSS) 
in  a  realistic  environment  near  an  active  trench  environment . An  earlier  test  of  the 
system  in  the  northwest  Pacific  had  been  unsuccessful  because  of  difficulties  in  drill¬ 
ing  an  acceptable  borehole  in  the  seafloor.  The  Scripps’  ship,  the  R/V  Melville,  was 
used  in  addition  to  the  D/V  Glomar  Challenger  for  tasks  related  to  site  surveying, 
ocean  bottom  seismograph  (OBS^  deployment  and  refraction  shooting.  Chief  Scientists 
on  the  expedition  were  John  Orcutt  and  Tom  Jordan  on  the  Melville  and  Bill  Menard  and 
Jim  Natland  on  the  Challenger.  The  experiment  was  quite  successful  and  accomplished 
all  major  goals.  MThe  MSS  was  recorded  aboard  the  D/V  Glomar  Challenger  for  several 
days  of  earthquake,  refraction  and  environmental  noise  experiments.  The  subsequent 
deployment  of  an  autonomous  recording  package  on  the  seafloor  was  successful,  -but  a 
leak  in  one  of  the  spheres  housing  the  batteries  resulted  in  a  power  loss  af ter"approx- 
imatcly  twn  days  The  ORS’s  launched  several  times  during  the  course  of  the  experiment 
operated  correctly  in  all  instances  including  the  45  day  teleseismic  recording  phase  at 
the  end.  The  MSS  recording  package  as  well  as  the  six  OBS’s  were  recovered  at  the  end 
of  the  experiment  by  the  R/V  Melville.  The  MSS  was  left  in  the  borehole  and  the  co¬ 
axial  cable  was  terminated  in  a  dummy  load.  The  mooring  system  used  for  recovery  was 
redeployed  in  the  event  a  future  experiment  became  necessary  or  desireable. 
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Chapter  1  of  this  report  discusses  a  propagation  model  for  the  oceanic  re¬ 
gional  phase  P  in  which  the  onset  consists  of  a  refraction  from  the  relatively 
high-Q  oceanic  lithosphere  while  the  extensive  coda  is  a  result  of  near-receiver 
sediment  and  water  reverberation.  This  model  contrasts  sharply  with  proposed 
scattering  mechanisms  for  generating  the  observed  coda  and  receives  strong  support 
from  the  numerical  modeling  and  the  high  quality  oceanic  earthquake  data  collected 
in  this  program.  In  addition  to  generating  realistic  coda  associated  with  P  and 
S  wavetrains,  the  model  has  successfully  predicted  the  absolute  amplitudes  of  ob¬ 
servations  collected  during  NGENDEI  so  that  the  model  can  be  confidently  used  for 
predictive  purposes.  The  reverberations  induced  in  the  sediment  and  ocean  column 
give  rise  to  "eigenf requencies"  associated  with  leaky  organ  pipe  modes  in  the  water 
column  (acoustic  only)  and  sediment  section  (elastic  waves).  These  eigenf requencies 
are  readily  isolated  in  both  the  synthetics  and  data  and  are  consistent  with  physical 
properties  measurements  conducted  at  the  drill  site.  The  horizontal  component  data 
and  synthetic  spectra  are  dominated  by  the  fundamental  shear  mode.  This  mode  gives 
rise  to  the  characteristic  monochromatic  reverberations  typical  of  marine  horizon¬ 
tal  components.  We  have  observed  that  the  P  phases  recorded  during  the  experiment 
are  quite  coherent  when  the  data  are  high  passed  above  the  fundamental  shear  eigen- 
frequency  for  sediment  reverberations.  This  is  a  sensible  behavior  given  that  the 
greatest  source  of  local  heterogeneity  on  the  seafloor  is  the  spatial  variation  of 
sediment  thickness.  Finally,  in  this  chapter,  we  note  that  the  observed  decrease 
in  P  group  velocity  at  a  range  of  about  30°  results  from  the  introduction  of  a 
mid-point  reflection  from  the  free  surface  of  the  lithosphere  refracted  phase  and 
that  the  required  thickness  of  the  lithosphere  is  quite  large  -  200  km. 

The  computations  which  support  the  work  on  the  regional  phases  are  quite  in¬ 
tensive  numerically  and  require  large  amounts  of  computer  time.  We  have  increased 
the  shear  Q  used  in  the  calculations  described  in  Chapter  1  from  225  to  1000  and 
have  repeated  the  synthesis  of  the  complete  seismograms.  The  S  phase  is  now  much 
more  realistic  than  that  presented  in  Chapter  1.  Computations  in  the  future  will 
be  run  with  the  goal  of  narrowing  the  current  loose  constraints  on  the  lithosphere 
attenuation  properties.  The  current  computation  is  being  extended  to  higher  fre¬ 
quencies  to  permit  a  more  thorough  study  of  the  Q.  The  algorithm  we  use  for  these 
complete  calculations  is  being  modified  to  run  on  a  CRAY  1  at  the  CRAY  headquarters 
in  Chippewa  Falls  and  will  be  used  extensively  on  the  CRAY  XMP/48  currently  being 
installed  at  the  NSF-sponsored  San  Diego  Supercomputer  Center. 

Chapter  2  presents  simultaneous  measurements  of  ocean  bottom  infrasonic  and 
ocean  bottom  and  sub-bottom  seismic  noise  in  the  frequency  band  0.1-20  Hz.  In 
the  band  0.1-1  Hz,  noise  propagation  as  seismic  modes  trapped  near  the  seafloor 
is  supported  by  observed  spectral  coherences,  cross-phases  and  ratios  between 
ocean  bottom  pressure  and  vertical  ground  motion,  and  by  the  relatively  lower 
noise  levels  in  the  borehole.  Noise  variations  in  this  band  are  clearly  cor¬ 
related  with  changes  in  local  wind  direction  and  speed,  presumably  through  ocean 
bottom  pressure  fluctuations  caused  by  nonlinear  wind  wave-wind  wave  and  wind  wave- 
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swell  interaction.  The  growth  of  spectral  peaks  and  generally  increased  levels 
over  the  microseism  band  in  both  seismic  and  acoustic  measurements  followed  sharp 
changes  in  wind  direction.  Activity  in  the  microseism  band  abated  as  the  wind 
persisted  or  died  away.  Amplitude  and  phase  relationships  between  pressure  and 
vertical  seismic  motion  at  the  microseism  periods  imply  that  a  portion  of  the 
noise  propagated  as  trapped  modes  of  the  ocean-seafloor  system.  The  appearance 
of  peaks  at  periods  corresponding  to  acoustic  reverberations  in  the  water-sediment 
column  suggests  that  some  of  the  ambient  noise  may  be  due  to  (leaky)  organ  pipe 
modes.  These  observations  constitute  the  first  comparison  of  simultaneous  ocean 
bottom  and  sub-bottom  seismic  noise  and  ocean  bottom  noise  over  the  microseism  band 
in  the  deep  ocean.  Recent  quantitative  experiments  with  the  same  algorithm  used  in 
the  wave  propagation  studies  discussed  above  have  resulted  in  a  successful  replication 
of  observed  noise  spectral  shapes  and  amplitudes. 

I  have  proposed  to  continue  aspects  of  this  research  upon  expiration  of  this 
grant  in  January,  1986.  In  particular,  I  have  suggested  comparing  noise  levels 
observed  on  the  seafloor  with  those  typically  encountered  on  islands.  Because  of 
breaking  waves,  islands  are  frequently  quite  noisy.  It  is  likely  that,  for  much 
of  the  discrimination  band,  adjacent  seafloor  noise  levels  are  much  lower.  If  this 
proved  to  be  the  case,  the  installation  of  stations  50-100  km  off  islands  with  co¬ 
axial  or  fiber  optic  cabling  could  enhance  current  detection  and  discrimination  ca¬ 
pabilities  . 

The  analyses  conducted  with  oceanic  regional  data  under  this  grant  have  con¬ 
centrated  on  the  results  of  the  NGENDEI  experiment.  I  have  proposed  to  check  ad¬ 
ditional  data,  in  particular  that  available  from  the  Wake  Island  Array,  for  the 
same  systematics  so  clearly  observed  in  the  NGENDEI  data.  This  study  would  serve 
as  a  rigorous  test  of  the  propagation  models  for  oceanic  regional  phases  which  we 
have  developed  during  the  past  year. 

Finally,  I  have  proposed  to  apply  many  of  the  things  we've  learned  about 
regional  wave  propagation  to  the  study  of  continental  phases  (L  ,  P  ,  R  ,  P  , 

S  )  as  recorded  by  the  new,  digital  NORESS  array.  These  phases  Ire,gof  §ourse, 
o?  great  interest  at  high  frequencies  in  nuclear  test  discrimination  and  defection. 

The  NORESS  array  is  ideally  suited  to  frequency  and  wavenumber  analyses  bit  such 
high  wavenumber  phenomena  and  comparisons  with  computed  responses.  The  oceanic 
data  are  intrinsically  limited  to  studies  in  the  temporal  and  frequency  domains 
since  the  small  number  and  large  inter-instrument  distances  preclude  wavenumber 
estimates.  I  feel  the  extension  of  these  quantitative  studies  to  the  NORESS  data 
would  be  immensely  rewarding  and  is  a  program  I  consider  to  be  quite  important. 


CHAPTER  I 


THE  PROPAGATION  OF  P. 


Thomas  Sereno  and  John  Orcutt 


Institute  of  Geophysics  and  Planetary  Physics 
Scripps  Institution  of  Oceanography,  A-025 
La  Jolla,  California  92093 
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ABSTRACT 

A  propagation  model  for  the  oceanic  PK  phase  is  presented  in  which  the  onset  con¬ 
sists  of  a  refraction  from  the  lower  oceanic  lithosphere  while  the  subsequent  coda  is  a 
result  of  near-receiver  sediment  and  water  reverberation.  The  first  arrival  of  the  oceanic 
Pn  phase  remains  entirely  within  the  lithosphere  and  for  epicentral  ranges  in  excess  of 
25  *  to  SO  *  contains  a  free  surface  reflection,  similar  to  the  seismic  phase  PP.  At  large 
epicentral  distances,  or  similarly  for  deep  events,  the  Pn  wavetrain  is  preceded  by  a  low 
frequency  forerunner  which  has  traversed  the  low  Q  asthenosphere.  The  coda  of  the  PH 
phase  was  analyzed  through  the  computation  of  synthetic  seismograms  using  wave 
number  integration.  Spectral  characteristics  associated  with  both  sediment  and  water 
reverberation  are  clearly  observed  in  both  the  synthetics  and  data  collected  in  the 
southwest  Pacific  during  the  1983  Ngendei  Seismic  Experiment. 

INTRODUCTION 

The  oceanic  P,  phase  is  characterized  by  low  spatial  attenuation,  a  high  frequency 
wavetrain  of  extremely  long  duration  and  an  onset  velocity  near  8.0  km/s1.  Through  the 
computation  of  synthetic  seismograms  in  a  laterally  homogeneous,  vertically  stratified 
medium  we  will  show  that  these  properties  can  be  achieved  in  the  absence  of  scattering, 
low  velocity  channels  or  high  Q  waveguides. 

SYNTHETIC  SEISMOGRAMS 

Wave  number  integration,  in  which  a  quadrature  over  wave  number  is  performed  at 
fixed  frequency  points,  was  used  to  generate  the  synthetic  seismograms3.  The  time 
domain  response  was  obtained  by  applying  a  Fourier  transform  to  the  resulting  complex 
spectrum.  The  seismograms  contain  all  reverberations,  including  those  off  the  free  surface, 
as  well  as  all  P  to  S  and  S  to  P  mode  conversions.  Attenuation  is  included  by  allowing 
the  material  velocities  to  be  complex  thereby  allowing  the  integration  to  be  performed 
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along  the  real  wave  number  axis.  The  most  serious  limitations  of  the  method  include  the 
restriction  to  models  consisting  of  a  stack  of  homogeneous  layers  and  the  high  computa¬ 
tional  costs  involved  for  models  with  more  than  10  to  15  layers. 

As  an  exploratory  model  we  have  adopted  that  of  Gettrust  and  Frazer*  which 
includes  a  4  km  thick  water  layer  and  a  500  m  thick  sediment  layer  overlying  a  vertically 
inhomogeneous  oceanic  crust  and  upper  mantle.  Compressional  and  shear  wave  velocities 
are  monotonically  increasing  functions  of  depth.  The  values  of  the  quality  factor,  Q,  were 
less  than  500  in  the  crust  and  upper  mantle  and  50,000  in  the  oceanic  water  column. 

Figure  1  is  an  example  of  complete  synthetic  seismograms  generated  for  a  sub-Moho 
thrust  fault  at  an  epicentral  range  of  1000  km.  The  late  arriving  energy  on  the  vertical 
component  (figure  la)  with  a  group  velocity  of  1.6  km/s  is  emergent  in  nature  and  has  a 
substantial  duration  characteristic  of  oceanic  abyssal  T-phases4.  Through  a  comparison  of 
delay  time  wavefields  between  Green's  functions  that  contain  this  late  arriving  energy 
and  Green’s  functions  that  do  not,  it  was  observed  that  the  ray  parameter  range  compris¬ 
ing  this  group  corresponds  to  that  of  post-critical  reflections  off  the  Moho.  The  synthetic 
T-phase  is  thus  explained  as  a  reverberation  between  the  free  surface  and  the  Moho  with 
a  substantial  number  of  water  and  crustal  reverberations  acting  to  extend  the  wavetrain. 
The  4.7  km/s  arrival  on  the  horizontal  component  (figure  lb)  is  grossly  identified  as  the 
oceanic  Sn  phase.  The  first  arrivals  of  this  group  are  horizontally  polarized  shear  waves 
and  are  thus  absent  from  the  vertical  component.  Work  is  continuing  in  an  effort  to  syn¬ 
thesize  more  realistic  5„  phases. 

The  synthetic  Pm  wavetrain  has  an  onset  velocity  of  8.1  km/s,  a  duration  of  over 
100  seconds  and  extends  into  the  onset  of  Sm .  The  amplitude  of  the  synthetic  PH 
wavetrain  was  within  a  factor  of  2  or  3  of  data  amplitudes  obtained  during  the  1983 
Ngendei  Seismic  Experiment  in  the  southwest  Pacific  (23*  49'S,  165*  32'W).  An  event 
located  near  the  Tonga  Islands  by  the  National  Earthquake  Information  Service  (NEIS) 
with  a  magnitude  of  5.8  and  a  moment  of  2.6  x  10*4  dyne-cm  was  recorded  by  three  ocean 
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bottom  seismometers  (OBS’s)  at  the  Ngendei  site.  The  maximum  amplitudes  of  these 
recordings  ranged  between  0.70  and  1.44  volts.  Convolution  of  the  synthetic  phase 
with  the  OBS  instrument  response  and  scaling  by  the  NEIS  published  moment  yielded  a 
maximum  amplitude  of  0.42  volts.  Differences  of  this  order  can  be  explained  simply  in 
terms  of  differences  in  model  parameters,  fault  orientation  and  focal  depth. 

REVERBERATIONS 

The  importance  of  reverberations  within  the  oceanic  water  column  and  sediment 
layer  is  most  effectively  demonstrated  in  the  frequency  domain.  It  is  simple  to  show  that 
the  power  spectrum  of  a  wavelet  reverberating  in  a  single  layer  has  eigenf’-equencies  given 
by 

/„=(  2n+l)-^-  n  =0,1,2,...  (1) 

where  dj  =  1/  T  and  T  is  the  time  between  multiples.  In  the  case  of  near-vertical 
incidence  T  can  be  approximated  by  the  two-way  travel  time  in  the  layer  and  (1) 
represents  the  eigenfrequencies  for  an  organ-pipe  mode*.  Arrivals  within  the  P„  wavetrain 
must  penetrate  the  mantle  in  order  to  have  group  velocities  between  6  and  8  km/s  and 
must  therefore  have  phase  velocities  greater  than  8  km/s.  A  result  of  this  is  that  the  rays 
are  quite  steep  at  the  receiver,  justifying  the  use  of  the  two-way  travel  time  in  equation 
(1).  The  «f//2  offset  of  the  spectral  peaks  is  due  to  a  * -phase  shift  between  multiples.  The 
width  of  the  spectral  peaks  at  frequencies  given  by  (1)  is  governed  by  the  impedance  con¬ 
trast  at  the  boundaries  of  the  layer.  Thus  leaky  modes,  which  suffer  a  continuous  loss  of 
energy  along  the  propagation  path  due  to  leakage  into  the  surrounding  medium,  will 
attenuate  more  rapidly  than  locked  modes  and  will  thus  have  substantially  wider  spectral 
peaks.  The  decay  of  the  reverberation  peaks  with  increasing  mode  number,  n,  is  a  func¬ 
tion  of  Q  in  the  layer. 


The  computation  of  the  full  synthetic  seismograms  pictured  in  figure  1  required 
extremely  fine  frequency  sampling  to  avoid  time  series  wrap-around.  The  calculation 
extended  to  a  Nyquist  frequency  of  only  6.4  Hz  which  involved  obtaining  the  response  at 
8193  frequency  points.  While  the  wave  number  integration  algorithm  is  not  strictly 
bandlimited,  the  costs  of  extending  the  calculation  to  higher  frequencies  are  formidable. 
Synthetic  power  spectra  do  not,  however,  require  this  fine  frequency  sampling  and  can  be 
computed  rather  efficiently  up  to  a  Nyquist  frequency  of  64.0  Hz.  Thus  direct  comparison 
between  data  collected  at  the  Ngendei  site  and  synthetics  is  easily  achieved  in  the  fre¬ 
quency  domain. 

The  Ngendei  site  was  the  location  of  hole  595B  of  Deep  Sea  Drilling  Project 
(DSDP),  Leg  91.  The  depth  of  the  oceanic  water  column  at  the  site  is  5.5  km  and,  from 
core  sample  measurements,  the  sediment  thickness  is  70  m  with  a  compressional  velocity 
of  1.6  km/s*.  For  compressional  wave  sediment  reverberation  the  two-way  travel  time  in 
the  layer  is  0.0875  s  yielding  an  eigenfrequency  spacing  of  df=  11  43  Hz.  The  impedance 
contrast  at  the  ocean-sediment  interface  is  quite  small  for  compressional  waves.  Thus  the 
sediment  layer  acts  as  an  extremely  leaky  waveguide  for  compressional  wave  propagation 
and  the  peaks  associated  with  sediment  reverberation  have  a  large  characteristic  width. 
The  reverberations  within  the  water  column  actually  include  the  sediment  layer  and  take 
place  between  the  free  surface  and  the  base  of  the  sediments.  The  resulting  eigenfre¬ 
quency  spacing  for  water  reverberations  is  therefore  d/=.  1375  Hz.  Both  interfaces  in  this 
case  have  very  strong  impedance  contrasts  and  the  associated  spectral  peaks  are  quite 
narrow.  Thus  the  power  spectrum  of  a  wavelet  reverberating  in  the  oceanic  water  column 
and  sediment  layer  consists  of  widely  spaced  spectral  "humps"  due  to  P-wave  sediment 
reverberation  superimposed  on  finely  spaced  spectral  "peaks"  due  to  water  reverberation. 

Figure  2  is  an  illustration  of  this  behavior.  Figure  2b  is  a  synthetic  power  spectrum 
computed  using  wave  number  integration  in  which  we  have  used  the  sediment  and  water 
properties  of  the  Ngendei  site.  In  an  effort  to  reduce  computational  costs  we  have 
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included  only  phase  velocities  between  6.9  and  25.0  km/s.  Figure  2a  depicts  the  first  3.0 
Hi  of  the  spectrum  with  vertical  lines  indicating  the  eigcnfrequencies  of  water  reverbera¬ 
tion  computed  using  equation  (1).  Figure  2b  contains  the  full  spectrum  after  correcting 
for  the  OBS  instrument  response  and  a  source  spectrum*  of  the  form 


(2) 


where  /  c  is  the  corner  frequency  and  was  given  a  value  of  0.6  Hz.  In  this  case  the  verti¬ 
cal  lines  indicate  positions  of  spectral  "humps"  associated  with  P-wave  sediment  rever¬ 
beration.  To  verify  our  interpretation  of  these  features  we  doubled  the  layer  thicknesses 
and  recomputed  the  spectra.  The  effect  was  to  decrease  the  spacing  of  the  spectral  peaks 
by  a  factor  of  2  thereby  validating  the  interpretation. 

Figures  2c  and  2d  contain  power  spectra  of  data  collected  at  the  Ngendei  site.  Fig¬ 
ure  2c  is  the  first  3.0  Hz  of  a  hydrophone  spectrum  illustrating  the  presence  of  water 
reverberations  and  figure  2d  is  an  OBS  vertical  component  spectrum  illustrating  the  pres¬ 
ence  of  sediment  reverberation.  These  spectral  features  were  quite  common  in  the  data 
and  were  independent  of  the  source  characteristics. 

The  presence  of  S-wave  sediment  reverberation  is  most  easily  detected  in  horizontal 
component  spectra  due  to  the  steep  incidence  of  the  rays.  For  shear  waves  the  impedance 
contrasts  at  both  the  ocean-sediment  and  sediment-crust  interfaces  are  quite  large  yield¬ 
ing  relatively  narrow  spectral  peaks.  The  low  quality  factor  for  shear  waves  in  the  sedi¬ 
ments,  however,  results  in  a  rapid  decay  of  the  spectral  peaks  with  increasing  mode 
number.  Figure  3a  depicts  the  first  3.0  Hz  of  the  synthetic  horizontal  component  power 
spectrum  after  including  the  OBS  instrument  response  and  source  spectrum  given  by  (2). 
The  shear  wave  sediment  velocity  used  in  the  computation  of  the  synthetic  spectrum  was 
116  m/s  yielding  an  eigenfrequency  spacing  of  0.83  Hz.  Figure  3b  is  the  synthetic  spec¬ 
trum  without  including  the  instrument  response.  Figure  Sc  is  the  first  3.0  Hz  of  the  power 
spectrum  of  horizontal  component  data  taken  at  the  Ngendei  site  and  Sd  was  obtained 
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after  removing  the  instrument  response.  In  figures  3a-3d,  dashed  vertical  lines  indicate 
predicted  frequencies  (equation  (l))  of  spectral  peaks  associated  with  S-wave  sediment 
reverberation.  The  positions,  widths  and  spectral  decay  rates  of  these  peaks  all  agree  with 
pre-dictions  for  shear  wave  sediment  reverberation.  Once  again  we  altered  the  sediment 
thickness  and  recomputed  the  spectra  to  verify  our  interpretation  of  the  peaks  in  the  syn¬ 
thetic  spectra. 


COHERENCE 

Strong  coherence  of  the  wavetrain  between  OBS’s  separated  by  large  distances 
is  not  expected  unless  the  water  and  sediment  properties  at  the  two  sites  are  nearly  equal. 
It  was  found  in  the  Ngendei  data  that  far  greater  coherence  could  be  achieved  by 
bandpassing  the  data  between  2.0  and  25.0  Hz  prior  to  comparison.  The  result  of  even  a 
small  change  in  either  S-wave  sediment  velocity  or  sediment  thickness  would  result  in  a 
measurable  shift  in  the  dominan .  frequency  of  the  wavetrain  (figure  3).  The  frequency 
band  between  0.0  and  2.0  Hz  is  the  most  sensitive  to  these  changes  resulting  in  a  reduc¬ 
tion  of  coherence  at  these  lower  frequencies.  Figure  4  is  a  comparison  of  data  collected  by 
two  different  OBS’s  separated  by  a  distance  of  1.9  km  for  the  same  event  after  bandpass¬ 
ing  between  2.0  and  25.0  Hz.  The  first  water  reverberation  is  evident  at  about  9  seconds 
and  the  excellent  coherence  well  into  the  P„  wavetrain  is  evidence  for  reverberation,  as 
opposed  to  scattering,  as  the  dominant  mechanism  of  P„  coda  generation.  Work  is  con¬ 
tinuing  on  the  problem  of  coherence  as  well  as  that  of  polarization  within  the  P„  coda. 

PROPAGATION  PATH 

The  observation  that  the  group  velocity  of  Pn  is  generally  higher  at  epicentral 
ranges  between  12*  and  30’  than  at  distances  outside  this  range7  can  be  explained  by 
the  confinement  of  P„  to  the  oceanic  lithosphere.  At  ranges  less  than  30  ’  the  onset  of  P„ 
is  associated  with  a  normal  refracted  P-wavc  remaining  in  the  lithosphere,  however  at  dis- 


tances  in  excess  of  SO  *  the  first  arrival  that  remains  in  the  lithosphere  contains  a  free  sur¬ 
face  reflection,  similar  to  the  seismic  phase  PP,  thereby  decreasing  the  group  velocity. 
This  is  demonstrated  in  figure  5b  which  displays  the  group  velocity  of  Pn  as  a  function  of 
epicentral  distance  for  the  velocity  model  shown  in  figure  5a.  In  order  to  achieve  a  cross¬ 
over  distance  of  30  * ,  a  lithosphere  thickness  of  nearly  200  km  is  required.  The  decrease 
in  group  velocity  at  ranges  less  than  12  *  is  the  result  of  shallower  turning  points  in  the 
lithosphere.  A  similar  interpretation  was  presented  by  Menke  and  Richards'  in  terms  of 
whispering  gallery  modes  at  the  Moho. 

Figure  5c  is  a  travel  time  curve  illustrating  the  arrival  times  of  normal  refracted  P 
verses  Pn  for  the  model  in  figure  5a.  The  curves  diverge  at  an  approximate  range  of  25  ‘ 
where  the  normal  refracted  arrival  is  a  ray  that  has  crossed  the  low  Q  asthenosphere  con¬ 
stituting  the  commonly  observed  low  frequency  forerunner  of  the  P„  wavetrain*.  Figure  6 
b  an  example  illustrating  the  normal  refracted  P  forerunner  of  oceanic  Pm .  For  deep 
events  the  separation  of  normal  refracted  P  and  Pm  occurs  for  smaller  ranges  as  PH  must 
propagate  up  and  around  the  corner  of  the  descending  slab.  During  the  Ngendei  Experi¬ 
ment  the  deployment  of  a  continuously  recording  borehole  seismometer  enabled  the 
recording  of  strong  Pm  phases  for  events  with  focal  depths  in  excess  of  450  km.  The 
suggestion  that  a  descending  slab  can  act  as  a  waveguide  to  propagate  the  high  frequency 
oceanic  Sn  phase  was  made  by  Isacks  and  Barasangi10.  Their  model  consbts  of  a  transfor¬ 
mation  of  the  basaltic  oceanic  crust  to  a  thin,  high  velocity,  eclogite  layer  as  the  slab  des¬ 
cends  to  greater  depths.  Thus  the  upper  flank  of  the  slab  provides  an  effective  reflecting 
surface  and,  in  combination  with  a  velocity  gradient  within  the  slab,  could  provide  an 
efficient  waveguide  for  the  propagation  of  high  frequency  sebmic  energy. 

SUMMARY 

The  model  of  the  propagation  of  Pn  presented  consists  of  an  arrival  confined  to  the 
oceanic  lithosphere  with  subsequent  water  and  sediment  reverberations  comprising  the 
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ensuing  wavetrain.  The  low  spatial  attenuation  observed  for  Pn  does  not  require  an  int¬ 
rinsically  high  Q  lithosphere  waveguide  but  rather  can  be  explained  in  terms  of  the 
extremely  high  Q  h  the  oceanic  water  column.  The  depletion  of  high  frequencies  in  the 
normal  refracted  arrival  relative  to  P„  is  due  to  propagation  through  the  low  Q  astheno- 
sphere. 
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Figure  1.  Complete  synthetic  seismograms  at  an  epicentral  range  of  1000  km  for  a 
thrust  fault  source  at  14.0  km  depth.  The  fault  had  a  dip  of  45’,  a  rake  or 
270  ’  and  an  azimuth  relative  to  the  receiver  of  00  * .  (a)  Vertical  component 
and  (b)  horizontal  component  rotated  25  *  clockwise  from  the  radial  direction. 
Group  velocities  of  8.1,  4.7,  2.6,  and  1.6  km/s  ore  indicated.  Phase  velocities 
between  S.85  and  25.0  km/s  were  included  in  the  calculation. 
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Figure  2.  Power  spectra  illustrating  the  presence  of  reverberation  within  the  Pn  coda. 

(a)  First  3.0  II*  of  the  synthetic  vertical  component  spectrum  depicted  in  (b). 
(c)  First  3.0  IIz  of  a  hydrophone  spectrum  and  (d)  a  complete  OBS  vertical 
component  spectrum.  In  figures  (a)  and  (c)  vertical  lines  indicate  predicted 
cigcnfrcqucncics  associated  with  water  reverberation  (equation  (l)).  In  figures 

(b)  and  (d)  vertical  linc3  indicate  predicted  cigenfrcquencics  for  P-wave  sedi¬ 


ment  reverberation. 
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Figure  3.  Power  spectra  illustrating  shear  wave  sediment  reverberations  within  the 
Pa  coda,  (a)  First  3.0  Hz  of  the  synthetic  horizontal  spectrum  after  including 
the  instrument  response  and  source  spectrum,  (b)  Synthetic  spectrum  of  (a) 
without  including  the  instrument  response.  Power  spectra  of  horizontal  com¬ 
ponent  data  taken  at  the  Ngcndci  site  (c)  before  removing  the  instrument 
response  and  (d)  after  removing  the  instrument  response.  In  (a)-(d),  dashed 
vertical  lines  indicate  predicted  eigcnfrequcncics  of  S-wave  sediment  reverbera¬ 
tion.  The  peak  amplitudes  in  each  figure  were  scaled  to  a  value  of  one. 
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Figure  5.  (a)  Velocity  model  used  to  compute  group  velocities  and  travel  times  of  Pn 
and  normally  refracted  P.  The  results  are  shown  for  a  surface  source  and 
receiver,  (b)  Group  velocity  of  P„  as  a  function  of  cpicentral  distance  for  the 
model  in  (a).  Group  velocity  has  been  defined  as  cpicentral  distance  divided 
by  total  travel  time.  The  discontinuous  break  at  SO  ’  represents  the  crossover 
between  a  ray  normally  refracted  in  the  lithosphere  and  one  that  contains  a 
free  surface  reflection,  (c)  Travel  time  curve  of  P„  (solid)  verses  normal 
refracted  P  with  turning  points  below  the  base  of  the  lithosphere  (dashed)  for 
the  model  in  (a). 
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Figure  6.  Vertical  component  seismogram  recorded  by  a  borehole  seismometer  at  the 
Ngendei  site.  The  focal  depth  of  the  event  was  293  km  in  the  Kermadec 
trench  (m4  =  5.0,  A=  12.91').  The  complex  P  waveform  between  2  and  4 
seconds  is  a  triplication  from  the  400  km  discontinuity.  Arrivals  following  5.5 
seconds  constitute  the  Pn  arrival.  The  enrichment  in  high  frequency  content 
-^n  relative  to  P  is  indicative  of  the  low  Q  of  the  sub-lithosphcrc  mantle. 
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ABSTRACT 

Simultaneous  measurements  of  ocean-bottom  infrasonic,  and  ocean-bottom 
and  sub-bottom  seismic  noise  in  the  frequency  band  0.1  -  20  Hz  are  presented. 
The  data  were  obtained  in  5.5  km-deep  water  in  the  South  Central  Pacific  with  a 
triaxial  borehole  seismograph  and  six  triaxial  ocean-bottom  seismographs  having 
externally-mounted  hydrophones.  The  borehole  sensors  were  emplaced  50  meters 
within  basement  rock  overlain  by  70  meters  of  pelagic  clay.  In  the  band  0.1-1 
Hz,  noise  propagation  as  seismic  modes  trapped  near  the  seafloor  is  supported  by 
observed  spectral  coherences,  cross-phases,  and  ratios  between  ocean-bottom  pres¬ 
sure  and  vertical  ground  motion,  and  by  the  relatively  lower  noise  levels  in  the 
borehole.  Noise  variations  in  this  band  are  clearly  correlated  with  changes  in  local 
wind  direction  and  speed,  presumably  through  ocean-bottom  pressure  fluctuations 
caused  by  nonlinear  wind  wave-wind  wave  and  wind  wave-swell  interaction. 


INTRODUCTION 

In  considering  the  response  of  the  coupled  ocean-seafloor  elastic  system,  one 
must  inescapably  conclude  that  ocean-bottom  seismic  and  infrasonic  noise  are 
related  in  source  and  character.  This  relationship  is  confirmed  by  observations  of 
these  two  types  of  ambient  noise.  Measurements  of  seismic  noise  at  the  ocean 
bottom  reveal  the  presence  of  the  ubiquitous  microseism  peak  (actually  an 
agglomeration  of  peaks)  ^over  the  microseism  band  (periods  between  1-10 
seconds)  (Bradner  et  al.  and  Latham  and  Sutton  ).  Ocean-bottom  pressure 
measurements  of  high  resolution  over  the  microseism  baud,  of  which  very  few  are 
published  (Latham  et  al.  ,  Nichols4,  and  Cox  et  al.  )  also  display  the  microseism 
peak.  Observed  amplitude  and  phase  relationships  between  the  ocean-bottom 
pressure  and  vertical-component  seismic  noise  are  consistent  with  those  of  normal 
modes  trapped  near  the  ocean-seafloor  interface  (Latham  and  Nowroozi  ). 

Over  the  microseism  band  an,  intimate  relationship  between  microseismic 
noise  and  non-linear  surfac^  gravity  wave  interactions  has  been  both 
hypothesized  (Longuet-IJjggins  and  Hasselmann  1  and  observed  (Haubrich  et 
al.  ,  Latham  and  Sutton  ,  Kibblewhite  and  Ewan*  and  Cox  et  al.  ).  This  non¬ 
linear  interaction  occurs  between  opposed  gravity  waves  of  the  same  period  and 
produces  ocean-bottom  pressures  fluctutions  of  half  the  gravity  wave  period.  The 
microseism  peaks  are  the  manifestation  of  this  phenomonenon. 

Although  regional  swell  is  most  frequently  identified  as  the  primary  source  of 
microseisms  in  the  deep  ocean,  recent  observations  have  emphasized  the  role 
played  by  local  wind-generated  waves  (Kibblewhite  and  Ewan1  ).  Hasselmann 
predicts  a  spectral  form  for  ocean-bottom  pressure  fluctuations  a  f  3  S2(f/  2)  at 
frequency  f,  where  S(f)  is  the  gravity  wave  spectrum.  The  asymptotic  form  of 
wi^-generated  waves  is  S(f)  oc  f~5  (Pierson  and  Moskowitz  *  ,  Hasselmann  et 
al.  )  or,  in  some  circumstances,  f-4  (Kahma  ).  jThus,  &  fall-off  oc  f~5—  f-7 
would  be  expected  and,  in  fact,  is  observed  (Cox  et  al.  ,  Latham  et  ah'*,  Figure  6, 
and  Nichols  j. 

This  paper  presents  ocean-bottom  and  sub-bottom  observations  from  a 
recent  ambient  noise  experiment  in  the  South  Central  Pacific  in  which  the  dom¬ 
inant  mechanism  of  ambient  noise  generation  was  clearly  wind-related.  The 


growth  of  spectral  peaks  and  generally  increased  levels  over  the  microseism  band 
in  both  seismic  and  acoustic  measurements  followed  sharp  changes  in  wind  direc¬ 
tion.  Activity  in  the  microseism  band  abated  as  the  wind  persisted  or  died  away. 
Amplitude  and  phase  relationships  between  pressure  and  vertical  seismic  motion 
at  the  microseism  periods  imply  that  the  a  portion  of  the  noise  propagated  as 
trapped  modes  of  the  ocean-seafloor  elastic  system.  The  appearance  of  peaks  at 
periods  corresponding  to  acoustic  reverberations  in  the  water-sediment  column 
suggests  that  some  of  the  ambient  noise  may  be  due  to  (leaky)  organ-pipe  modes. 

These  observations  constitute  the  first  comparison  of  simultaneous  ocean- 
bottom  and  sub-bottom  seismic  noise  and  ocean-bottom  acoustic  noise  over  the 
microseism  band  in  the  deep  ocean. 
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I.  EXPERIMENT  DESCRIPTION 

The  Defense  Advanced  Research  Projects  Agency  (DARPA)  developed  a 
marine  borehole  seismograph,  the  Marine  Seismic  System  (MSS),  for  seismic 
discrimination  applications.  The  successful  test  of  an  MSS  prototype  tethered  to 
the  D/V  Glomar  Challenger  on  Deep  Sea  Drilling  Project  (DSDP)  Leg  78B  (Bal¬ 
lard  et  al.^)  led  to  deployment  of  the  fully-operational  MSS  on  DSDP  Leg  91, 
which  composed  the  first  portion  of  the  Ngendei  Experiment. 

A  primary  objective  of  the  Ngendei  Experiment  was  the  appraisal  of  the 
relative  seismic  signal  and  noise  characteristics  at  a  marine  borehole  site  and 
proximate  ocean-bottom  site.  This  behavior  was  at  issue  because  the  first  gen¬ 
eration  of  ocean-bottom  seismograph  (OBS)  efforts  during  the  1960’s  revealed 
that  ambient  seismic  noise  levels  in  the  microseism  band  were  typically  10-20  dB 
higher  at  ocean-bottom  sites  than  at  nearby  land  sites.  In  the  microseism  band, 
seismic  noise  at  land  sites  is  dominated  by  fundamental  mode  surface  waves,  the 
amplitude^,  of  which  diminishes  with  depth  into  more  rigid  basement  rock 
(Sorrells1  ).  Since  marine  seismic  noise  is  similarly  composed,  it  was  thought  that 
the  same  advantage  might  obtain  at  a  marine  borehole  site  over  an  ocean-bottom 
site. 

The  MSS  sensors  were  deployed  in  DSDP  Hole  595B,  located  midway 
between  New  Zealand  and  Tahiti  at  28'  49.34'  S,165  *  31.61'  W.  Further 
details  of  the  Ngendei  Experiment  may  be  found  in  Volume  88/91  of  the  Initial 
Reports  of  the  DSDP. 

n.  INSTRUMENTATION  AND  OPERATIONS 


A.  Ocean  Bottom  Seismographs 

The  six  OBS’s  are  modifications  of  the  fou ^-component,  digital  event- 
recorders  described  by  Moore  et  al.^  (Orcutt  et  al.^  ).  Output  from  three  com¬ 
ponents  of  ground  motion  (one  vertical  and  two  horizontal)  and  a  hydrophone 
are  amplified,  bandpassed,  and  digitized  with  a  12-bit  word  at  128 
samplcs/second  (s/s).  The  gain  of  each  channel  is  2N,  where  N=0,l,...,9.  A 
channel’s  gain  exponent  N  is  modified  only  when  the  long-term  average  of  its 
digital  output,  corrected  for  DC  offset,  falls  outside  pre-defmed  bounds  of  a  few 
least  significant  bits. 


Semiconductor  memory,  nominally  13-seconds  long  per  channel,  continuously 
accumulates  the  four  channel  data  in  a  cyclically  overwritten  buffer  to  avoid  the 
mechanical  shaking  and  electrical  transients  of  recorder  operation.  Software 
writes  events  to  tape  at  scheduled  times  or  when  the  short  term  average  of  verti¬ 
cal  component  amplitudes  exceeds  a  long  term  average  threshold.  Additional 
RAM  boards  extended  the  13-second  buffers  of  several  OBS’s  to  59  seconds. 


The  OBS’s  employed  two  seismometer  responses  (see  Figure  la),  selectable 
only  prior  to  deployment.  The  ‘teleseismic’  response  emphasizes  lower  frequen¬ 
cies  more  than  the  ‘refraction’  response  in  order  to  increase  the  sensitivity  of  the 
OBS  event-triggering  algorithm  to  moderate-amplitude,  high-frequency  depleted 
earthquake  signals.  The  ‘refraction’  and  ‘teleseismic’  capsules  employed  the  same 
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hydrophone  response  ( Figure  la). 

The  seismometers  are  the  1-Hz  free  period  Mark  Products  L-4-3D  triaxial 
set,  and  the  hydrophone  is  an  Ocean  and  Atmospheric  Sciences,  Inc.,  model  E- 
2PD  deep-water  sensor.  A  bouyant  aluminum  sphere  mounted  on  a  detachable, 
dish-footed  tripod  houses  the  OBS  instrumentation.  Explosive  bolts,  detonated 
either  by  an  acoustic  or  programmed  command,  free  the  instrumentation  sphere, 
which  then  floats  to  the  surface  for  recovery. 

B.  Marine  Seismic  System 

The  MSS  comprises  a  borehole  instrumentation  package  (BIP)  connected  by 
coaxial  cable  to  either  an  ocean-bottom  recording  unit,  known  as  the  Bottom 
Processing  Package  (BPP),  or  to  two  shipboard  recording  systems  employed  for 
diagnostic  checkout  and  recording  redundancy.  One  of  the  shipboard  recording 
systems,  referred  to  as  the  Teledyne  system,  specifically  monitors  the  borehole 
package,  while  the  other,  the  Gould  system,  emulates  the  BPP.  A  special  tool  at 
the  end  of  the  drill  string  carries  the  sensor  package  to  the  seafloor,  although  sim¬ 
ple  lowering  of  the  mechanically-strong  coaxial  cable,  ’wireline’  deployment,  pro¬ 
vides  an  alternate  means  to  achieve  this  end.  To  date,  only  drill  string  deploy¬ 
ments  have  been  employed,  although  successful  recovery  of  the  BIP  has  been 
achieved  by  wireline. 

The  borehole  sensors  are  four  Teledyne  Model  S-750’s  configured  as  a  pri¬ 
mary  triple  (orthogonal  horizontals,  X  and  Y,  and  one  vertical  Z)  and  a  backup 
vertical  component  (B).  The  backup  sensor  B  is  located  175.5  cm  below  the  pri¬ 
mary  vertical  sensor  Z.  Output  from  the  three  primary  sensors  X,  Y,  and  Z  are 
filtered  in  mid-period  (MP)  and  short-period  (SP)  bands,  while  output  from  the 
backup  sensor  B  is  passed  through  the  SP  band  only.  The  MP  and  SP  transfer 
functions  are  shown  in  Figure  lb.  The  resultant  seven  data  streams  are  then 
digitized  and  transmitted  to  recording  devices  via  the  coaxial  cable. 

Data  recorded  on  the  Teledyne  system  are  continuously  gain-ranged  to  yield 
a  21-bit  dynamic  range  with  14-bit  resolution.  No  low-pass  filter  is  employed 
other  than  the  roll-off  of  the  instrument  response.  The  MP  and  SP  bands  are 
digitized  at  4  and  40  s/s,  respectively,  although  either  of  the  two  SP  vertical 
channels  may  be  digitized  at  80  s/s  to  the  exclusion  of  the  other. 

Output  recorded  by  the  BPP  and  its  shipboard  counterpart  are  oversampled 
256-fold,  transmitted,  digitally-convolved  with  a  low-pass  filter,  resampled  at  the 
proper  rates  of  4  s/s  and  40  s/s  for  the  MP  and  SP  bands,  respectively,  and  buf¬ 
fered  to  tape  as  24-bit  integers  with  24-bit  resolution.  The  transfer  function  of 
the  digital  low-pass  filter  is  nearly  flat  to  75%  of  the  Nyquist  frequency,  and 
thereafter  falls  at  approximately  720  dB/decade.  The  low-pass  filter  coefficients 
are  the  same  for  both  passbands.  The  dashed  lines  in  Figure  lb  illustrate  the 
augmentation  due  to  the  low-pass  filter 
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C.  Operations 

The  first  20  days  of  the  60-day  long  Ngendei  Experiment  comprised  a 
seismic  refraction  experiment,  during  which  the  OBS’s  recorded  refraction  signals 
and  teleseismic  signals  with  two  distinct  arrays,  and  the  MSS  operated  with  its 
two  shipboard  systems.  The  remaining  40  days  were  devoted  to  a  teleseismic 
experiment  during  which  the  OBS’s  formed  a  single  array  for  earthquake  and 
noise  recording,  and  the  MSS  employed  its  ocean-bottom  recording  system. 

The  D/V  Glomar  Challenger  conducted  MSS  installation  and  shipboard 
recording  operations.  The  Scripps  Institution  of  Oceanography  R/V  Melville 
tended  the  six  OBS’s  and  recovered  all  instrumentation.  Figures  2  and  3  show 
instrumentation  locations  throughout  the  experiment.  In  these  figures,  octogons 
centered  on  ellipses  depict  estimated  OBS  locations  and  their  95%  confidence 
bounds  based  on  Ieas^quares  modeling  of  satellite  and  acoustic  ranging  data 
(Creager  and  Dorman1  ).  Diamonds  and  arrows  denote  OBS  drop  points  com¬ 
puted  from  satellite  navagation  and  Melville  radar  ranging  of  Challenger.  The 
borehole  location  is  computed  relative  to  the  best  satellite  fix  of  Challenger’s 
location. 

During  the  refraction  experiment,  four  OBS’s  were  placed  0.6-0.7  km  from 
the  MSS  site  ( Figure  2a),  and  two  at  approximately  30  km  (Figure  2b).  The 
former  quartet  of  capsules  employed  the  ‘refraction’  seismic  response,  13-second 
event  memory,  and  scheduled  event  recording,  while  the  latter  pair  employed  the 
‘teleseismic’  response,  59-second  event  memory,  and  triggered  event  recording 
with  scheduled  noise  samples  at  three-hour  intervals. 

The  BIP  rested  unclamped  at  the  bottom  of  DSDP  Hole  595B,  approxi¬ 
mately  50  meters  within  a  basaltic  basement  overlain  by  a  70-meter  thick  layer  of 
zeolitic  pelagic  clays  bearing  hard  stringers  of  porcellanite  and  chert  near  its  base. 
The  BIP  holelock/cable  isolation  mechanism  damaged  the  coaxial  cable  during  an 
aborted  deployment  and  was  subseqeuntly  discarded.  Comparison  of  MSS  and 
OBS  waveforms  indicate  that  MSS  sensor  coupling  with  the  ground  was  ade¬ 
quate.  A  coil  of  the  coaxial  cable  payed  out  on  the  seafloor  aided  in  cable-motion 
isolation. 

During  the  teleseismic  experiment,  the  six  OBS’s  were  redeployed  in  a  tri¬ 
angular  array  with  single  capsules  forming  two  25-km  legs  at  approximate  right- 
angles  to  the  south  and  east  of  the  MSS  site  (Figure  3a).  The  other  four  capsules 
were  clustered  within  3  km  to  the  west  and  north  of  the  MSS  in  a  trapeziform 
array  (Figure  Sb).  All  OBS’s  employed  the  ‘teleseismic’  response,  and  triggered 
event  recording,  two  with  13-second  memory  buffers  and  four  with  59-second 
buffers.  Daily  noise  events  were  recorded  at  Universal  Time  (UT)  midnight.  The 
MSS  ocean-bottom  recording  system  failed  two  days  after  entering  the  water  due 
to  a  water  leak. 

Successful  recovery  of  all  OBS’s  and  the  MSS  recording  system  took  place 
from  Melville  over  the  period  22-26  March,  1983.  Following  a  two-day  shipboard 
recording  period,  the  MSS  recovery  mooring  was  redeployed  with  a  dummy  load 
to  replace  the  recording  system.  Both  BIP  and  recovery  mooring  are  still  in  place 
for  future  use. 
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HI.  NOISE  DATA 

Totals  of  119,  43,  and  83  hours  of  seismic  data  were  obtained  from  the 
Teledyne  system,  the  Gould  ocean-bottom  recording  system,  and  its  shipboard 
counterpart,  respectively.  During  the  refraction  experiment,  114  and  70  hours  of 
shipboard  data  were  recorded  on  the  Teledyne  and  Gould  systems,  respectively. 
OBS’s  recorded  a  total  of  1400  noise  events  at  740  distinct  times,  140  con¬ 
currently  with  MSS  data. 

Various  MSS  data  channels  were  unavailable  at  certain  times.  The  horizon¬ 
tal  sensor  Y  malfunctioned  throughout  the  entire  Ngendei  experiment  for  unk¬ 
nown  reasons  and  provided  no  usable  data.  A  shipboard  source  of  electrical 
interference  contaminated  the  subsequently-disconnected  Gould  channel  SZ.  A 
Gould  electronics  failure  in  its  BIP  bay  disabled  Gould  channel  MZ  at  the  start 
of  the  teleseismic  experiment.  Finally,  as  mentioned  above,  one  of  the  Teledyne 
SP  Z  channels  was  not  available  whenever  the  other  recorded  at  80  s/s. 


IV.  DATA  QUALITY 

Data  quality  was  assessed  with  both  time-  and  frequency-domain  criteria, 
the  latter  applied  to  power,  cross,  and  coherence  spectra. 

Rj|iable  spectral  esthetes  were  computed  in  the  manner  described  by 
Welch  and  Carter  et  a.1.  .  Non-overlapping,  de-meaned  segments  of  digital 

data  were  tapered  with  a  normalized  Hanning  window  and  Fourier  transformed. 
The  resulting  spectra  were  averaged  to  achieve  a  reduction  of  estimate  variance. 
Power  spectra,  defined  only  for  non-negative  frequencies,  are  normalized  so  that 
the  total  power  in  the  Nyquist-frequency  bandwidth  equals  the  time-domain  vari¬ 
ance,  a  normalization  which  essentially  folds  the  power  of  negative  frequencies 
into  the  positive. 

Estimates  of  statistical  reliability  were  calculated  under  the  assumption  that 
the  data  derived  from  a  stationary,  Gaussian  stochastic  process.  Power  spectral 
estimates  are  then  distributed  as  x2  random  variables  with  2N  degrees  of  free¬ 
dom,  where  N  is  tfye^ number  of  spectra  averaged  to  obtain  the  estimate  (see,  e.g. 
Bendat  and  PiersoF1). 

The  complicated  statistics  of  the  magnitude-squared  coherence,  here  termed 
coherence  for  conciseness^^are  essentially  those  of  the  linear-regression  correlation 
coefficient  (Carter  et  al.  ;  Fischer  ).  The  coherence’s  cumulative  probability 
function  simplifies  greatly  if  the  actual  coherence  C  is  zero,  from  which  the  coher¬ 
ence  R  at  which  an  estimate  C  is  significantly  different  from  zero  for  a  given  sig¬ 
nificance  level  is  easily  determined.  This  is  inferred  from 

P(C  <  R  |  C  =  0  )  =  1  -  (1  -  R)n-\  (1) 

where  P(— )  is  the  probability  that  the  estimate  C  does  not  exceed  the  value  R 
given  the  the  actual  coherence  C  and  N  is  the  number  of  spectra  averaged  in  the 
coherence  estimate.  Thus,  the  threshold  coherence  R  is  given  at  the  100P  %  con¬ 
fidence  level  by 


R  =  1  -  (1  -  P)  N“ 


(2) 
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In  computing  spectral  estimates  from  OBS  data,  the  number  of  spectra  aver¬ 
aged,  and  stationarity  must  trade-off  because  OBS  noise  samples  are  available 
only  at  discrete  times:  N  must  be  large  enough  to  yield  estimate  stability  and 
useful  error  bounds,  but  not  so  large  as  to  sample  over  a  time  interval  such  that 
the  noise  field  has  statistically  evolved.  N=10,  a  30-90  minute  interval,  was  found 
to  be  generally  optimal.  This  trade-off  is  not  very  important  for  MSS  data  since 
they  were  continuously  recorded  except  during  brief  interuptions.  However,  where 
comparisons  of  simultaneous  observations  are  desired,  the  statistics  and  frequency 
resolution  of  the  MSS  estimates  were  made  similar  to  those  of  the  OBS. 

Linear  regressions  of  Gould  data  against  Teledyne  data  in  conjunction  with 
sensor  calibrations  revealed  a  3  dB  shortfall  of  Gould  amplitudes.  The  calibra¬ 
tions  also  show  that  instrument  response  below  0.2  Hz  is  not  good  for  any  MSS 
channel. 

Figure  4  compares  power  spectral  density  estimates  of  simultaneous  vertical- 
component  Teledyne  and  Gould  data  with  instrument  noise  spectra  which  have 
not  been  corrected  for  instrument  response.  The  peaks  above  8  Hz  in  the  SB 
spectra  are  ship-generated.  Peaks  in  the  SB  Teledyne  spectrum  (Figure  4C) 
which  are  absent  in  the  SB  Gould  spectrum  (Figure  4b)  are  manifestations  of 
aliasing  in  the  Teledyne  data  due  to  inadequate  high-frequency  filtering.  Gould 
and  Teledyne  spectra  have  comparable  values  near  the  Nyquist  frequency,  indi¬ 
cating  probable  malfunctioning  of  Gould’s  low-pass  filter. 

Also  shown  in  Figure  4  are  estimates  of  system  noise,  provided  by  H.  B. 
Durham  of  Sandia  Laboratories.  The  system  noise  is  not  significant  over  the 
entire  SP  passband  for  either  Teledyne  or  Gould  data.  The  flat  spectrum  of 
quantizing  noise  for  the  Teledyne  data,  associated  with  uniformly -^distributed 
round-off  errors  of  digital  representation  (see  e.g.  Bendat  and  Piersoi  pp.  231- 
232),  has  a  level  near  10  L  counts  /Hz,  well  below  the  Teledyne  signal  level. 
Estimates  of  Gould’s  digitizing  noise  spectrum,  that  of  a  delta  modulator,  are  not 
available.  The  MZ  spectrum  (Figure  4a)  is  similarly  untainted  by  system  and 
sampling  noise.  The  MZ  instrument  response  rolls-off  sufficiently  at  high  frequen¬ 
cies  to  avoid  aliasing. 

Spectral  coherences  ( Figure  5)  further  delineate  frequencies  usable  for 
ambient  noise  estimates.  Coherences  exceeding  0.28,  illustrated  in  Figure  5,  with 
a  dash  line,  are  considered  significant  at  the  95  %  level.  MP  coherences  (Figure 
5a)  typically  exceed  0.95  below  the  low-pass  filter  knee  at  1.5  Hz  except  at  1  Hz, 
where  it  drops  sharply  over  a  very  narrow  band  of  frequencies.  The  sharp  drop  is 
due  to  energy  present  in  the  Gould  data  but  not  the  Teledyne,  and  may  be 
related  to  the  Gould  digitizing  scheme.  Gould  and  Teledyne  SB  coherences  (Fig¬ 
ure  5b)  exceeds  0.8  between  approximately  0.2  and  4.0  Hz.  The  coherence  at 
higher  frequencies  declines  amidst  peaks  due  to  ship-generated  signals.  The  sharp 
drops  at  approximately  2  and  4  Hz  are  similar  to  that  at  1  Hz  in  the  MP  coher¬ 
ence.  The  frequencies  bounding  the  SP  band  of  high  coherence  varied  a  great 
deal  during  the  MSS  deployment,  but  generally  were  those  shown  in  Figure  5. 
Coherence  estimates  between  SB  and  SZ  Teledyne  data  (Figure  5c)  were  similar. 
Differences  at  frequencies  greater  than  4  Hz  are  probably  due  to  the  combined 
effects  of  aliasing  and  actual  signal  differences  over  the  175.5-cm  sensor 
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separation.  The  high  coherence  between  the  Gould  and  Teledyne  data  at  the 
Nyquist  frequencies  indicates  a  possible  deficiancy  in  the  Gould  low-pass  filtering 
scheme. 

OBS’s  employing  the  refraction  response  can  measure  lower-amplitude 
ground  motion  at  high  frequencies  than  those  employing  the  teleseismic  response. 
The  refraction  response  was  specifically  designed  to  ‘whiten’  the  ambient  signal  in 
the  presence  of  typical  deep-sea  microseismic  noise  (Moore  et  al.A  ),  while  the 
teleseismic  response  was  designed  to  low-pass  it.  Filtered  high-frequency  ampli¬ 
tudes  are  consequently  less  than  that  represented  by  a  least  significant  bit  as 
recorded  with  the  teleseismic  response.  Figure  6  compares  spectra  measured  with 
the  different  OBS  responses  (upper  curves)  and  shows  the  corresponding  spectra 
of  least-count  noise  (lower  curves).  The  spectra,  corrected  for  instrument 
response,  were  not  measured  simultaneously.  Although  most  of  the  data  compos¬ 
ing  the  teleseismic  capsule’s  spectrum  (solid  upper  curve)  were  recorded  after  ship 
departure,  their  spectral  levels  exceed  ship  noise  levels  measured  with  a  refraction 
capsule  (dashed  lower  curve)  at  most  frequencies.  In  the  teleseismic  spectrum, 
sampling  noise  clearly  dominates  at  frequencies  greater  than  15  Hz,  and  is  prob¬ 
ably  significant  in  the  range  7-18  Hz. 

The  lowest  valid  frequency  of  OBS  microseismic  measurements  is  taken  to  be 
0.1  Hz.  At  lower  frequencies,  least-count  and  system  noise  dominate.  The  system 
noise  was  inferred  from  the  lack  of  estimate  stability  over  a  wide  range  of 
numbers  of  spectra  averaged  to  compute  power  density  estimates. 

Ship-generated  noise  must  be  estimated  on  a  case-by-case  basis  because  it 
depends  on  ship  speed,  activity  and  range.  Generally,  sharply  distinct  peaks 
superimposed  on  a  broadband  level  characterize  ship-generated  noise.  In  the 
Ngendei  data,  most  peak  frequencies  appear  to  be  multiples  of  som^  fundamental, 
indicating  that  propellor  cavitation  caused  the  peaks  (e.g.  Uricx  ).  Challenger 
generated  peaks  at  multiples  of  10  Hz  while  maintaining  station,  and  two  others 
at  approximately  6.5  and  7.5  Hz  while  under  way.  Melville  generated  peaks  at 
multiples  of  6  and  9  Hz.  Figure  7  shows  the  growth  of  a  peak  at  18  Hz  due  to 
Melville’s  proximity  to  the  MSS. 

Figure  8a  illustrates  ship-generated  noise  with  a  comparison  of  vertical- 
component  power  density  spectra,  corrected  for  gain  but  not  instrument  response, 
of  OBS  Phred  data  recorded  before  (solid  curve)  and  after  (dashed  curve)  both 
ships  departed.  The  capsule,  located  1.6  km  west  of  the  borehole,  employed  the 
teleseismic  response.  The  measurement  prior  to  ship  departure  was  made  as  both 
Challenger  and  Melville  hove  to  within  2  km  of  the  OBS.  In  this  instance,  ship¬ 
generated  noise  clearly  dominates  the  spectrum  at  all  frequencies  above  approxi¬ 
mately  4.5  Hz.  After  ship  departure,  least-count  noise  dominated  above  7  Hz, 
although  the  lowest  frequency  at  which  it  dominated  varied  by  ±  1  Hz  since  the 
level  depends  on  OBS  gain. 

Some  ambiguity  exists  concerning  ship  noise  in  the  MSS  data.  Because  the 
BPP  failed  two  days  after  its  deployment  when  Challenger  was  13  km  away  dril¬ 
ling  Hole  596  and  Melville  was  en  route  to  Tahiti,  MSS  noise  measurements  in 
the  complete  absence  of  ships  are  not  available.  However,  Challenger  was  briefly 
silent  c  0510Z,  13  Feb,  while  monitoring  an  ocean-  bottom  navigational 


-  29- 


transponder,  thus  providing  a  sample  of  low  ship  noise  for  the  MSS.  Figure  8b 
compares  the  MSS  spectrum  of  this  time  period  with  that  measured  a  short  time 
later  when  Challenger  began  using  its  dynamic  positioning  system  thrusters.  Also 
shown  in  Figure  8b  is  an  estimate  of  SP  system  noise  which  has  been  scaled  to  fit 
the  ‘quiet’  spectrum.  The  good  fit  suggests  that  above  7  Hz  system  noise  dom¬ 
inates  the  spectrum. 

The  dominance  of  ship-generated  noise  is  in  many  instances  diagnosed  by  an 
abrupt  change  in  spectral  slope,  such  as  is  seen  in  Figure  8  at  approximately  4 
Hz.  On  this  basis,  it  is  apparent  from  these  and  other  spectra  that  ship  noise  may 
be  important  at  frequencies  as  low  as  2  or  3  Hz  in  both  the  OBS  and  MSS  data. 

In  summary,  reliable  spectral  estimates  of  ambient  seismic  noise  are  avail¬ 
able  from  MSS  MP  data  between  0.2  and  1.5  Hz;  from  MSS  SP  data  between  0.2 
and  7  Hz;  and  from  OBS  data  between  0.1  and  7  Hz.  These  bands  are  subject  to 
the  variability  described  above  due  to  ship  proximity,  and  OBS  gain  range  and 
instrument  response.  At  frequencies  greater  than  7  Hz,  sampling  noise  dominates 
OBS  teleseismic  data,  so  that  lower  bounds  on  ocean-bottom  microseismic  noise 
levels  must  be  established  using  data  recorded  with  the  refraction  response. 

V.  RESULTS  AND  DISCUSSION 

Figures  9-12  compare  displacement  power  densities  for  all  available  com¬ 
ponents,  in  units  of  nanometers2/  Hz  (nm2/  Hz)  for  simultaneous  MSS  and  OBS 
data.  The  data  were  obtained  at  5-  and  10-minute  intervals  over  a  70-minute 
period.  Spectra  are  averages  of  10  individual  spectra  obtained  from  approximately 
13-seconds  of  data.  The  OBS,  Suzy,  was  0.6  km  north  of  the  MSS  site.  Based  on 
considerations  mentioned  above,  ship  noise  is  assumed  to  dominate  these  spectra 
above  4  Hz. 

The  spectra  all  display  the  ubiquitous  microseism  peak  near  0.2  Hz,  where 
the  observed  noise  levels  are  comparable,  approximately  108  nm2/  Hz.  At  higher 
frequencies,  ambient  ocean-bottom  microseismic  noise  levels  exceed  the  borehole 
levels  by  5-30  dB  ( Figures  9  and  10).  Common  to  the  spectra,  though  not  prom¬ 
inent  in  all  of  the  OBS  spectra,  are  peaks  centered  on  0.5  and  0.7  Hz.  At  the 
ocean  bottom,  horizontal  levels  exceed  the  vertical  levels  at  all  frequencies  greater 
than  that  of  the  microseism  peak  and  display  a  character  quite  different  from 
that  of  the  vertical  component  between  0.3  and  0.G  Hz  (Figure  11).  In  contrast, 
the  borehole  vertical-component  spectrum  uniformly  exceeds  the  horizontal- 
component  spectrum  by  an  insignificant  amount  (4  dB)  over  the  entire  frequency 
band  (Figure  12).  These  observations  are  consistent  with  the  propagation  of 
noise  as  a  fundamental  mode  interface  wave  with  large  vertical  and  horizontal 
displacements  at  the  interface  and  lower  amplitudes  at  depth  in  more  rigid  rock. 
The  comparable  levels  at  the  microseism  peak  suggest  that  the  amplitude  of  the 
trapped  waves  are  still  significant  at  the  depth  of  the  borehole  sensor. 

The  ocean-bottom  pressure  spectrum  measured  simultaneously  with  the 
above  spectra  also  displays  the  microseism  peak  (Figure  IS).  The  fall-off  towards 
higher  frequencies  decreases  with  frequency  from  the  microscism  peak.  Between 
0.2  and  0.7  it  falls  at  100  dB/decade;  between  0.7  and  2  Hz  at  54  dB/dccade;  and 
above  2  Hz  at  27  dB/decade,  where  ship  noise  evidently  dominates. 
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The  microseism  peak  and  a  0.4  Hz-wide  band  centered  at  1  Hz  display  high 
coherence  and  a  cross-phase  of  approximately  180'  (Figure  1\  (a)  and  (b)) 
between  the  pressure  and  vertical  displacement.  It  may  be  shown  that  at  a 
fluid-solid  interface,  kinematic  boundary  conditions  require  a  180 '  phase  shift  for 
both  no^ally-incident  plane  waves  ^reverberating  in  the  water  column 
(Bradner  ,  and  Schneider  and  Backus^  )  and  for  trapped  modes  of  the  ocean- 
seafloor  elastic  system,  the  gravest  o^  which  is  confined  near  the  interface  at  these 
frequencies  (Latham  and  Nowroozi  ,  and  Latham  and  Sutton;.  For  plane 
waves,  the  pressure-seismic  velocity  ratio  is  constant  with  frequency  and  equals 
the  product  of  density  p  and  acoustic  speed  a  in  the  water,  1.5xl06  in  MKS 
units,  while  for  normal  modes  it  may  vary  with  frequency  and  exceed  pa  (Lat¬ 
ham  et  al.  and  Latham  and  Nowroozi  ).  Figure  14  (c)  plots  the  ratio  between 
the  pressure  and  vertical-component  velocity  power  density  spectra,  normalized 
by  (/ja)2.  Over  the  two  bands  of  high  coherence,  this  ratio  is  on  the  order  of  10, 
rather  than  one,  suggesting  that  in  this  instance  they  represent  trapped  seismic 
modes.  Thus,  in  this  case  the  practice  of  comparing  seismic  ground  velocities 
V(f)  to  a  ‘plane-wave  equivalent’  pressure  P(f)  via  the  relation  P=/>cV  is 
incorrect. 

Figures  15-16  shows  the  coherence  and  cross-phase  between  the  vertical  and 
hydrophone  channels  of  OBSs  Suzy  and  Lynn,  seperated  by  0.6  km,  for  the  same 
time  interval  as  above.  Again,  energy  at  the  microseism  peak  is  highly  coherent. 
The  cross-phase  at  these  frequencies  is  not  significantly  different  from  zero, 
implying  either  simultaneous  arrival  of  microseism  energy  at  the  two  capsules  or 
a  wavelength  much  larger  that  the  capsule  seperation.  Latham  et  al.6  calculate 
the  dispersion  for  normal  modes  in  an  ocean-bottom  structure  similar  to  that  of 
the  Ngendei  site.  At  these  frequencies  the  fundamental  mode  has  a  phase  velo¬ 
city  Cp  of  1.0  -  1.5  km/s,  corresponding  to  a  wavelength  A=Cp/  f  on  the  order  of 
5  km.  The  capsule  separation,  roughly  A/  10,  is  too  large  to  account  for  the 
observed  small  cross-phases,  which  rules-out  oblique  arrival  of  the  microseism 
energy  from  a  distant  source.  It  therefore  appears  that  the  microseism  source 
location  in  this  instance  is  either  local  or  along  a  line  that  is  approximately  nor¬ 
mal  to  a  line  defined  by  the  OBSs. 

Evidence  abounds  that  the  sea  state  affects  microscismic  and  infrasonic  noise 
levels.  The  non-linear  interaction  of  opposed  ocean  surface  waves  of  a  given  fre¬ 
quency  generate  substantial  ocean^bottom  pressure  fluctuations  of  twice  the 
waves’  frequency  (Longuet-Higgins  and  Hasselmann^).  This  has  been  con¬ 
clusively  demonstrated  for  swell  reflected  from  a  shoreline  (Haubrich  et  al.^),  and 
for  wind-waves  (Kibblcwhite  and  Ewan  ;  Webb  ).  Dccklog  observations  pro¬ 
vide  observations  pertinent  to  local  and  distant  noise  source  discrimination.  Fig¬ 
ure  19  summarizes  Challenger  decklog  observations  of  wind  and  swell  between  24 
January  and  13  February,  1983  (Julian  dates  24-44)  at  two-hour  intervals  except 
for  8-10  hours  swell  gaps  corresponding  to  local  nighttime,  and  for  hourly  obser¬ 
vations  during  days  33-37.  The  prevalent  swell  arrived  from  the  south  and  had 
6-7  second  periods.  The  highest-amplitude  swell  had  periods  between  8-10 
seconds,  predominantly  9  seconds,  and  arrived  from  the  south  between  dates  29- 
31  and  34-35.  Very  little  northerly  swell  was  observed.  The  wind  direction 
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changed  both  suddenly  and  gradually,  the  latter  being  marked  by  smooth 
changes  in  wind  speed. 

Figure  18  shows  the  variation  of  OBS  Phred  vertical-component  acceleration 
spectra  with  time  at  three-hour  intervals  except  for  gaps  due  to  triggered  event 
recording.  The  spectra,  single  realizations  corrected  for  instrument  response,  are 
computed  from  58  seconds  of  data.  The  noise  clearly  changes  with  time,  particu¬ 
larly  at  the  prominent  peaks  near  0.2  and  0.5  Hz,  which  have  comparable  levels. 
Figure  18(b)  plots  the  time-variation  of  the  power  in  a  0.2  Hz  band  centered  on 
these  peaks,  along  with  wind  ( Figure  19(a) )  and  swell  ( Figure  19(c) )  observations 
over  the  same  period.  Abrupt  increases  in  the  power  of  both  peaks,  e.g.  days  28.0, 
30.0  and  33.0,  coincide  in  several  instances  with  sudden  changes  in  wind  direc¬ 
tion.  The  onset  of  high-amplitude  9-second  swell  coincides  with  the  wind  change 
and  power  increases  at  day  30.0,  but  opposing  9-second  swell,  probably  masked 
by  other  swell,  were  not  observed.  The  onset  at  day  33  of  the  high-speed  wind 
episode  from  the  SSE  coincides  with  a  large  power  transient  of  the  0.5  Hz  peak 
and  a  gradual  ramping  of  power  of  the  0.2  Hz  peak.  Noise  observations  of  13- 
seconds  length  collected  at  three-minute  intervals  from  the  ‘refraction’  OBS’s 
near  the  MSS  site  during  day  33  show  the  peak  at  0.5  Hz  building  quickly  over 
the  course  of  an  hour.  Power  for  both  peaks  gradually  diminished  as  the  wind 
weakened  and  wound  northward.  OBS  Phred  pressure  spectra  varied  similarly 
with  time.  The  clearly  correlated  behavior  of  the  ambient  noise  and  the  wind  as 
a  major  source  of  microseisms. 

The  observed  variations  of  the  wind  are  clearly  related  to  those  of  the  noise 
in  Figure  19.  A  sharp  change  in  the  wind  vector  immediately  preceeds,  to  within 
the  2-hour  resolution  of  the  dccklog,  noi|e  increases  at  frequencies  greater  than 
the  microseism  peak  band.  Hasselmann  shows  that  the  direct  conveyance  of 
pressure  from  the  atmosphere  at  these  frequencies  is  unlikely,  leaving  only  the 
non-linear  interaction  of  wind-generated  gravity  waves  or  of  wind-waves  with 
background  swell  as  a  likely  source^for  the  seismic  variations.  In  the  case  of  wind 
wave  interaction  only,  Hasselmann  derives  the  following  expression  for  the  spec¬ 
tral  density  Sp(f}  of  the  ocean-bottom  pressure  fluctuations  at  frequency  f: 

2  2 

Sp(2f)  =  8jt  w3S^r(f)  f G [9  )G(0  -Hr)  d9  ,  (3) 


where  w  =  27rf,  Sw(f)  is  the  wave  height  power  density,  p  and  c  are  the  density 
and  sound  speed  of  water,  g  is  gravitational  acceleration  at  sea  level,  and  G(0  )  is 
the  angular  dependence  of  the  wind  wave  spectrum.  Sp(f)  may  be  thought  of  as  a 
forcing  function  for  seismic  wave  generation.  Note  that  the  non-linear  interaction 
between  exactly-opposed  waves  of  frequency  f  (cf  the  angular  integrand  in  Equa¬ 
tion  (S))  leads  to  the  double-frequency  dependence  2f  of  the  ocean-bottom  pres¬ 
sure  fluctuations.  Implicit  in  Equation  (8)  is  a  white  wavenumber  spectrum 
which  contributes  only  at  phase  speeds  exceeding  the  sound  speed  of  water,  that 
is  for  wavenumbers  ^u;/  c. 

jgA  large  body  of  observations  (Pierson  and  Moskowitz**;  Hasselmann  et 
al.  )  support  seif-similar  form  for  the  wind  wave  spectrum  Sw(f)  of 
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Sw(f)  a  r5  e  4  1  , 


(4) 


The  function  described  in  Equation  (4)  falls  exponentially  for  f^fm  and  as  f~5  for 
f»fm,  where  fm  is  the  peak  frequency.  Empirically,  fm  decreases  and  power  den¬ 
sity  increases  with  time  under  the  influence  of  a  constant  wind.  Tyler  et  al. 
found  that  G(0  )  is  adequately  described  by 

cos  q(0  /  2),  (5) 

where  q  is  a  function  of  frequency  and  0=0  corresponds  to  the  mean  wind  direc¬ 
tion. 


A  more  detailed  view  of  this  interplay  may  be  obtained  from  the 
continuously-collected  MSS  data.  Figures  20  (a),  (b)  show  the  temporal  variation 
of  mid-period  vertical  acceleration  in  the  borehole  using  data  from  Teledyne 
channels  MZ.  Averages  of  14  spectra  computed  from  consecutive  128-second 
samples,  nearly  30  minutes  total  data,  form  the  individual  spectra  of  Figure  20. 
The  dominant  feature  of  these  spectra  is  the  microseism  peak  near  0.2  IIz,  which 
with  time  steadily  migrated  towards  lower  frequencies,  and  grew  larger  in  ampli¬ 
tude  and  narrower  in  width.  Note  the  ‘saddle’  in  the  peak’s  amplitudes  near  day 
41  coinciding  with  with  generally- increased  levels  at  higher  frequencies.  In  con¬ 
trast  to  the  ocean-bottom  accelerations,  levels  at  the  0.5  Hz  peak  ( Figure  18  (b)) 
are  an  order  of  magnitude  lower  than  those  at  0.2  Hz.  Other  peaks  are  observed 
at  0.35,  0.65,  and  0.80  Hz  ( Figure  20(b)). 

Figure  21(b)  compares  the  temporal  variation  of  power  for  the  microseism 
(‘primary’)  and  0.5  Hz  ancillary  peaks  of  Figure  20.  Power  was  estimated  by 
integrating  over  0.10  Hz-wide  bands  centered  on  0.22  and  0.50  Hz.  Power  varia¬ 
tions  in  these  two  bands  are  roughly  parallel  until  c.  day  40.7,  after  which  they 
were  complimentary,  corresponding  to  the  microseism  peak’s  ‘saddle’.  The  power 
of  these  peaks  decreased  in  an  unknown  manner  early  in  the  MSS  shipboard 
recording  period,  then  began  growing  c.  day  39.0.  Power  in  the  ancillary  band 
diminished  during  day  39,  but  resurged  c.  day  40.0.  Power  in  both  bands  changed 
little  in  the  half-day  preceding  day  40.7,  then  decreased  in  the  primary  band  and 
increased  in  the  ancillary  band.  At  c.  day  41.0,  power  leveled-off  in  the  microse¬ 
ism  peak  band  and  generally  increased  at  higher  frequencies,  which  peaked  c.  day 
41.2.  The  rate  of  power  growth  in  the  two  bands  increased  c.  day  41.5,  so  that 
primary  power  increased  and  ancillary  power  fall-off  slowed.  By  the  end  of  the 
shipboard  recording  period,  power  in  both  bands  attained  levels  comparable  to 
those  observed  c.  day  40.5. 

Figures  21  (a)  and  (c)  show  weather  observations  during  the  MSS  shipboard 
data  collection  period,  which  began  just  as  the  epoch  of  the  strongest  winds 
observed  during  the  Ngendei  Experiment  was  concluding.  Southerly  swell  of  6-7 
second  period  dominated  most  of  the  time  although  northerwesterly  swell  briefly 
dominated  c.  day  41  (Figure  21  (e)).  The  wind,  on  the  other  hand,  varied  a 
great  deal  ( Figure  21  (a)).  A  moderate  northeasterly  wind  blew  during  most  of 
days  37  and  38  at  a  nominal  10  knots.  During  day  39,  the  wind  was  confused, 
though  still  northeasterly.  The  wind  blew  northwesterly  throughout  day  40,  and 
southwesterly  for  the  rest  of  the  shipboard  recording  period.  The  highest  wind 
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speeds  during  the  MSS  shipboard  recording  period,  18  knots,  accompanied  the 
shift  in  wind  direction  from  northwest  to  southwest  at  day  41.0. 

The  available  data  permit  only  a  gross  description  of  the  sea.  However,  the 
forgoing  discussion  of  wind-wave  evolution  permits  drawing  reasonable  inferences 
of  the  wind-wave  behavior  from  Challenger  decklog  observations.  The  increases  in 
noise  power  at  frequencies  greater  than  those  of  the  microseism  band  which  fol¬ 
low  soon  after  sharp  changes  in  the  wind  vector  are  probably  due  to  the  wind 
degrading  the  pre-existing  wave  field.  The  subsequent  power  increases  at  these 
frequencies  evident  in  Figure  21,  c.  days  39.0,  40.0,  and  41.0  are  due  to  the 
newly-evolving  wind  wave  field. 

The  behavior  of  the  microseism  peak  band  is  not  as  easily  explained  by  the 
available  data.  The  growth  of  power  throughout  days  38  and  39,  as  well  as  the 
shift  towards  lower  frequencies  apparent  in  Figure  20(a)  suggests  an  evolving 
wind  wave  source.  Similarly,  the  fall  during  day  40  and  subsequent  rise  during 
day  41  suggests  the  demise  of  the  old  wave  regime  and  domination  of  the  new 
regime  in  response  to  winds  which  blew  approximately  90 0  ,  then  180  *  with 
respect  to  the  old  wind. 


If  wind  waves  were  soley  responsible  for  the  observed  seismic  noise  varia¬ 
tions,  Equation  (4)  implies  that  a  fall-off  of  the  wind-wave  spectrum  would  result 
in  an  ocean-bottom  pressure  spectrum  <x  {~7,  rather  than  the  observed  f~5  [Fig¬ 
ure  IS).  However,  the  ragged  fall-off  observed  in  all  spectra  suggests  that  the 
seafloor  responded  resonantly  to  the  wind-wave  forcing.  Sharp  peaks  at  approxi¬ 
mately  0.35,  0.5,  0.65,  and  0.8  Hz  were  observed  in  the  borehole  data.  It  will  be 
recalled  that  the  peaks  at  0.5  and  0.65  Hz  were  prominent  in  one  or  more  of  the 
OBS  spectra  in  Figure  11.  Equal  spacing  of  spectral  peaks  is  expected  for  so- 
called  ‘organ  pipe’  modes,  due  to  the  reverberation  of  acoustic  waves  between  the 
sea  surface  and  the  sediment-basement  interface.  Such  reverberations  are  an 
important  constituent  of  the  seismic  phase  Pn  observed  in  the  ocean  (Sereno  and 
Orcutt  )  which  in  this  context  is  comprised  of  P  waves  travelling  almost 
entirely  in  the  upper  mantle  before  turning  sharply  to  the  vertical  in  the  sedi¬ 
ment.  The  gravest  peaks  are  at  frequencies  (2n+l)  df/2,  n  =  0,  1,  ...,  and  are 
therefore  seperated  by  df,  where  df  is  the  two-way  travel  time  of  a  single  rever¬ 
beration.  At  the  MSS  site^  ,  df  =  0.14,  very  near  the  observed  0.15  Hz.  The 
nonlinear  interaction  of  gravity  waves  described  above  doqbtless  generates  acous¬ 
tic  waves^in  the  cmnpressibl^  ocean  (Longuet-Higgins  ,  Hasselmann8,  Bre- 
khovskikh  ,  Hughes  ,  Lloyd0  ),  so  that  a  significant  contribution  to  the 
observed  ambient  noise,  both  acoustic  and  seismic,  may  be  due  to  these  rever¬ 
berations.  On  the  other  hand,  if  the  ocean-bottom  pressure  fluctuations  due  to 
the  wind  waves  are  spatially  restricted,  they  must  perforce  possess  energy  over  a 
broad  range  of  wavenumbers,  and  therefore  be  capable  of  exciting  trapped  seismic 
modes  which  have  a  much  smaller  phase  velocity  than  the  "leaky"  acoustic  rever¬ 
berations.  Modelling  the  relative  contributions  to  the  ambient  noise  of  these 
complementary  wave  types  will  be  the  topic  of  further  research. 


VI.  CONCLUSIONS 


The  observations  reported  here  constrain  the  source  mechanisms  and  propa¬ 
gation  modes  of  infrasonic  and  microseismic  noise  in  the  deep  ocean.  High  coher¬ 
ences,  180  *  cross-phases,  and  ratios  between  displacement  and  pressure  at  fre¬ 
quencies  below  1  Hz  imply  a  trapped  wave  propagation  mode.  Modelling  of  these 
measurements,  as  well  as  the  horizontal-component  measurements,  may  further 
constrain  the  wave  type,  although  a  definitive  experiment  to  measure  the 
wavenumber  spectrum  may  be  required. 

Episodic  variations  of  the  wind  direction  and  speed  correlate  with  changes  in 
both  acoustic  and  microseismic  noise.  In  particular,  a  peak  at  0.5  Hz  responded 
immediately  to  sudden  changes  in  wind  direction,  both  in  ocean-bottom  and  sub¬ 
bottom  observations.  The  microseisin  peak  is  also  affected  by  the  wind,  an  obser¬ 
vation  that  is  in  contrast  to  multitudes  of  previous  reports  focussing  attention  on 
swell  as  the  chief  agent  of  noise  generation.  The  mechanism  for  transferring 
energy  from  wind  to  microseisms  is  presumably  the  non-linear  interaction  of 
wind-generated  waves.  Spectral  shapes,  and  temporal  evolution  of  pressure  and 
displacement  spectra  support  this  hypothesis. 
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Displacement  responses  of  (a)  ocean-bottom  seismograph  (OBS)  and  (b) 
Marine  Seismic  System  (MSS)  instrumentation  used  during  the  Ngendei 
experiment.  (a)  OBS  seismometer  and  hydrophone  responses.  Two 
responses  were  available  for  the  OBS’s  (ocean  bottom  seismographs),  (b) 
The  Could  MSS  mid-period  (MP)  and  short-period  (SP)  responses  differ 
from  t lie  Teledyne  responses  by  the  low-pass  FIR  (Finite  Impulse  Response) 
digital  filter  shown  with  dashed  lines. 


Locations  of  the.  MSS  and  six  OBS’s  during  the  refraction  experiment.  See 
text  for  explanation  of  symbols,  (a)  Configuration  of  the  four  capsules  which 
recorded  refraction  signals  (Juan,  Karen,  Lynn,  and  Suzy).  (b)  Configuration 
of  the  two  OBS’s  which  concurrently  recorded  teleseismic  signals  in  a  trig¬ 
gered  mode  (Janice  a.nd  Phred). 


(a)  and  (b):  Locations  of  the  OBS’s  during  the  teleseismic  experiment.  See 
text  for  explanation  of  symbols. 


Spectral  power  density  estimates  of  typical  MSS  noise  samples  compared  to 
system  noise  levels.  The  spectra,  given  in  counts  /Hz,  have  not  been 
corrected  for  instrument  response,  (a)  Gould  Channel  MZ.  (b)  Gould  chan¬ 
nel  SB.  (c)  Teledyne  channel  SB. 


Coherence  spectra,  corrected  for  alignment,  between  (a)  Teledyne  and 
Gould  channel  MZ,  (b)  Teledyne  and  Gould  channel  SB,  and  (c)  Teledyne 
channels  SZ  and  SB. 


Comparison  of  noise  spectrum  obtained  with  the  teleseismic  (solid  upper 
curve)  and  refraction  (dashed  upper  curve)  OBS  responses.  The  correspond¬ 
ing  instrument-corrected  sampling  noise  spectra  (lower  curves)  are  also 
shown. 


Figure  7  Demonstration  of  Melville's  ship  noise  during  the  period  034:0430-035:0000 
(Julian  date),  (a)  Smoothed,  integrated  OBS  Suzy  vertical  component  power 
density  between  17.9  and  18.2  Hz.  {I)) Melville's  range  from  the  MSS  site. 


Figure  8  Illustration  of  ship  noise  in  the  (a)  OBS  and  (b)  MSS  data.  Solid  spectra 
were  observed  when  ships  were  nearby  and  operational,  and  the  dashed  spec¬ 
tra  when  ships  v/ere  completely  absent  (OBS)  or  still  nearby  (MSS).  MSS 
data  are  not  available  at  any  time  when  ships  were  completely  absent. 


Figure  9  Simultaneoxis  vertical-component  displacement  power  densities  at  the  ocean 
bottom  (dashed  line)  and  in  the  borehole  (solid  line). 


Figure  10  Simultaneous  horizontal-component  displacement  power  densities  at  the 
ocean  bottom  (dashed  line)  and  in  the  borehole  (solid  line). 

Figure  11  A  comparison  of  simultaneous  vertical  (solid)  and  horizontal  (dashed)  power 
densities  in  the  borehole. 


Figure  12  A  comparison  of  simultaneous  vertical  (solid)  and  horizontal  (dashed)  power 
densities  at  the  ocean-bottom. 


Figure  13  Ocean-bottom  pressure  measurements  concurrent  with  the  displacement 
spectra  of  Figures  9-12. 


Figure  14  (a)  Coherence  and  (b)  cross-phase  between  OBS  Suzy  channels  Z  (vertical) 
and  P  (hydrophone).  The  cross-phase,  plotted  in  multiples  of  7r,  has  been 
corrected  for  instrument  response  phase,  (c)  The  ratio  between  OBS  Suzy 
pressure  and  vertical  ground  velocity  power  densities,  normalized  by  the  pro¬ 
duct  of  density  and  acoustic  speed  in  the  water. 


Figure  15  (a)  Coherence  and  (b)  cross-phase  between  the  vertical  channels  of  OBS’s 
Suzy  and  Lynn,  as  for  Figure  14. 


Figure  1G  (a)  Coherence  and  (b)  cross-phase  between  the  hydrophone  channels  of 
OBS’s  Suzy  and  Lynn,  as  for  Figure  14. 


Figure  17  Challenger  decklog  observations  of  (a)  wind  and  (b)  swell  depicted  as  vectors 
with  origin  at  the  time  of  measurement,  lengths  proportional  to  wind  speed 
(knots)  or  swell  height  (feet),  and  direction  that  of  the  measurement.  The 
swell  measurements  are  plotted  at  a  level  corresponding  to  swell  period. 


Figure  18  Variation  with  time  of  OBS  Phred  vertical-component  acceleration  power  in 
the  frequency  band  0.1  -  1.0  Hz. 


Figure  19  Temporal  variation  of  (b)  peaks  at  0.2  and  0.5  Hz  of  the  OBS  Phred  data 
shown  in  Figure  18  compared  with  (a)  wind  and  (c)  swell  observations  shown 
in  Figure  17. 


Figure  20  Temporal  variation  of  vertical  accleration  in  the  borehole  in  the  frequency 
bands  (a)  0.1  -  1.0  Hz  and  (b)  0.4  -  0.9  Hz. 


Figure  21  Temporal  variation  of  (b)  peaks  at  0.2  and  0.5  Hz  of  the  MSS  data  shown  in 
Figure  20  compared  with  (a)  wind  and  (c)  swell  observations  shown  in  Fig¬ 
ure  17. 
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